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ABSTRACT 

We determine the distances to 18 galaxy clusters with redshifts ranging from z ~ 0.14 to z ~ 0.78 
from a maximum likelihood joint analysis of 30 GHz interferometric Sunyaev-Zel'dovich effect (SZE) 
£N) ' and X-ray observations. We model the intracluster medium (ICM) using a spherical isothermal f3 model. 

We quantify the statistical and systematic uncertainties inherent to these direct distance measurements, 
and we determine constraints on the Hubble parameter for three different cosmologies. These distances 
CNI ' imply a Hubble constant of 6OI4 t\g km s _1 Mpc -1 for an CIm = 0.3, J1a = 0.7 cosmology, where the 

£>Y uncertainties correspond to statistical followed by systematic at 68% confidence. With a sample of 18 

clusters, systematic uncertainties clearly dominate. The systcmatics are observationally approachable 
and will be addressed in the coming years through the current generation of X-ray satellites (Chandra 
& XMM- Newton) and radio observatories (OVRO, BIMA, & VLA). Analysis of high redshift clusters 
detected in future SZE and X-ray surveys will allow a determination of the geometry of the universe 
CNI , from SZE determined distances. 

Subject headings: cosmic microwave background — cosmology: observations - distance scale - galaxies: 
clusters techniques: interferometric 



1. INTRODUCTION 
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CO ■ Analysis of Sunyaev-Zel'dovich effect (SZE) and X-ray data provides a method of directly determining distances to 
I galaxy clusters at any redshift. Clusters of galaxies contain hot (k B T e ~ 10 keV) gas, known as the intracluster medium 
. (ICM), trapped in their potential wells. Cosmic microwave background (CMB) photons passing through a massive cluster 
' interact with the energetic ICM electrons with a r ~ 0.01 probability. This inverse Compton scattering preferentially 
boosts the energy of a scattered CMB photon causing a small (< 1 mK) distortion in the CMB spectrum, known as the 
• Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1970, 1972). The SZE appears as a decrement for frequencies < 218 
Q-T GHz and as an increment for frequencies > 218 GHz. The SZE is proportional to the pressure integrated along the line 
of sight AT cx / n e T e d£. X-ray emission from the ICM has a different dependence on the density S x cx J n^,A eH d£, where 
£^ ■ A eH is the X-ray cooling function. Taking advantage of the different density dependencies and with some assumptions 
^ \ about the geometry of the cluster, the distance to the cluster may be determined. SZE and X-ray determined distances 
C$ • are independent of the extragalactic distance ladder and provide distances to high redshift galaxy clusters. The promise 
' of direct distances has been one of the primary motivations for SZE observations. 

In the last decade, SZE detections have become routine due to advances in both instrumentation and observational 
strategy. Recent high signal-to-noise ratio detections have been made with single dish observations at radio wavelengths 
5h ' (Birkinshaw & Hughes 1994; Hcrbig ct al. 1995; Myers et al. 1997; Hughes & Birkinshaw 1998), millimeter wavelengths 
(Holzapfel et al. 1997a, b; Pointecouteau et al. 1999, 2001) and submillimetcr wavelengths (Lamarre et al. 1998; Komatsu 
et al. 1999). Interferometric observations at centimeter wavelengths are now routinely producing high quality images of 
the SZE (Jones et al. 1993; Grainge et al. 1993; Carlstrom et al. 1996; Carlstrom et al. 1998; Saunders et al. 1999; Grainge 
et al. 2002b; Reese et al. 2000; Grego et al. 2000; Carlstrom et al. 2000; Grego et al. 2001; Jones et al. 2002). 

In this paper, we present a maximum likelihood joint analysis of our 30 GHz interferometric SZE observations with 
archival Rontgen Satellite (ROSAT) X-ray imaging observations. Cluster X-ray temperatures, metallicity, and H I column 
densities are taken from the literature. The intracluster medium (ICM) is modeled as a spherical isothermal /3 model. 
We refine the analysis technique described in Reese et al. (2000) and apply it to a sample of 18 clusters for which we 
determine distances. These distances are then used to measure the Hubble constant. This is the largest homogeneously 
analyzed sample of SZE clusters with distance determinations thus far. To date, there are about 20 published estimates 
of Hq based on combining X-ray and SZE data for individual clusters (see Birkinshaw 1999, for a review and compiled 
distances). Most notably, those results include one sample consisting of 7 nearby (z < 0.1) galaxy clusters (Mason et al. 
2001; Mason 1999; Myers et al. 1997) and a sample of 5 intermediate redshift (0.14 < z < 0.3) clusters (Jones et al. 2002). 

The cluster sample selection for this paper is discussed in §2. The centimeter wave SZE system and interferometric 
SZE data are described in §3.1. A brief overview the ROSAT X-ray cluster data is given in §3.2. The analysis method, 
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including uncertainty estimation, is outlined in §4 along with the model fitting results. Distances and our determination 
of the Hubble parameter appear in §5. Sources of possible systematic uncertainties are discussed in §6. Section 7 contains 
a discussion of the results and future prospects. Throughout this work, the galaxy cluster CI 0016 + 16 (z = 0.546) will 
be used as the example cluster to illustrate both the analysis method and general results. All uncertainties are 68.3% 
confidence unless explicitly stated otherwise. 

2. CLUSTER SAMPLE 

The determination the Hubble parameter requires a large sample of galaxy clusters free of selection effects. For 
example, clusters selected by X-ray surface brightness will preferentially include clusters elongated along the line of sight. 
A spherical analysis will underestimate the line of sight length of the cluster, causing the derived Hubble parameter to 
be biased low. In theory, selecting by X-ray luminosity, Lx, alleviates the selection bias problem. In practice, this is 
complicated by the fact that X-ray surveys are surface brightness limited; clusters just at the detection limit that arc 
elongated along the line of sight will be detected while clusters just at the detection limit that are instead extended in the 
plane of the sky will be missed. Staying well above the detection limit of the survey will alleviate this potential pitfall. 
Observational considerations for cluster sample selection include the declinations of the clusters, the size (redshift) of the 
cluster, possible radio point sources in the cluster field, and SZE brightness for which we use as an indicator. The 
Owens Valley Radio Observatory (OVRO) and Berkeley-Illinois-Maryland Association (BIMA) interferometers have been 
used to observe known X-ray clusters with z > 0.14, declination > —15°, and L x > 5 x 10 44 h§ erg s _1 (0.3-3.5 keV 
band for Einstein and 0.1-2.4 keV band for ROSAT). In addition, short, preliminary observations of many clusters are 
also performed to investigate possible point sources in the field. 

The OVRO/BIMA SZE imaging project initially chose targets from the limited number of known X-ray bright clusters. 
With the publishing of X-ray cluster surveys, the OVRO/BIMA SZE imaging project chose targets from three X-ray 
catalogs of galaxy clusters: the Einstein Observatory Extended Medium Sensitivity Survey, EMSS (Gioia et al. 1990; 
Stocke et al. 1991; Gioia & Luppino 1994; Maccacaro et al. 1994); the ROSAT X-ray Brightest Abell Clusters, XBACs 
(Ebeling et al. 1996b,a); and the ROSAT Brightest Cluster Sample, BCS (Ebeling et al. 2000a; Crawford et al. 1999; 
Ebeling et al. 1998, 1997). We have also recently included two more recent R OSA T samples of distant massive clusters 
to our cluster selection database: the Wide Angle ROSAT Pointed Survey, WARPS (Fairley et al. 2000; Ebeling et al. 
2000b; Jones et al. 1998; Scharf et al. 1997); and the MAssive Cluster Survey, MACS (Ebeling et al. 2001). So far, we 
have high S/N detections in 21 clusters with redshifts z > 0.45. 

The distance calculation requires three data sets: SZE, X-ray imaging, and X-ray spectroscopic data. We have obtained 
high signal-to-noise detections of the SZE in 45 galaxy clusters. The subsample of these clusters that also have high 
signal-to-noisc X-ray imaging data and published electron temperatures contains 18 galaxy clusters. Table 1 summarizes 
the redshifts and X-ray luminosities for each galaxy cluster in our sample. 

3. DATA 

Here we briefly describe the SZE and X-ray observations and data reduction. Table 2 summarizes the observation 
times for both the SZE and ROSAT observations of each cluster in our sample. The SZE observation times are the total 
on-source integration times for the interferometric SZE data used in this analysis. The ROSAT observation times are the 
total livetimes of the pointings used in this analysis. 

3.1. Interferometric SZE Data 

The extremely low systematics of interferometers make them well suited to study the weak SZE signal. A unique 
feature of interferometers is their ability to separate the diffuse, negative SZE emission from small scale, positive point 
source emission through the spatial filtering of the interferometer. Interferometers also provide a well defined angular and 
spectral filter, which is important in the analysis of the SZE data discussed in §4. 

3.1.1. Centimeter- Wave System and Observing Strategy 

Over the past several summers, we outfitted the Berkeley-Illinois-Maryland Association (BIMA) millimeter array in 
Hat Creek, California, and the Owens Valley Radio Observatory (OVRO) millimeter array in Big Pine, California, with 
centimeter wavelength receivers. Our receivers use cooled (~ 10 K) High Electron Mobility Transistor (HEMT) amplifiers 
(Pospicszalski et al. 1995) operating over 26-36 GHz with characteristic receiver temperatures of T rx ^11-20 K over the 
28-30 GHz band used for the observations presented here. When combined with the BIMA or OVRO systems, these 
receivers obtain typical system temperatures scaled to above the atmosphere of T sys ~ 35-45 K. Most telescopes are 
placed in a compact configuration to maximize sensitivity on angular scales subtended by distant clusters (~ 1'), but 
telescopes are always placed at longer baselines for simultaneous detection of point sources. Every half hour we observe 
a bright quasar, commonly called a phase calibrator, for about two minutes to monitor the system phase and gain. The 
total integration time for each cluster field is given in Table 2 for both OVRO and BIMA. 

An interferometer samples the Fourier transform of the sky brightness distribution multiplied by the primary beam 
rather than the direct image of the sky. The SZE data files include the positions in the Fourier domain, which depend 
on the arrangement of the telescopes in the array and the declination of the source, the real and imaginary Fourier 
components, and a measure of the noise in the real and imaginary components. The Fourier conjugate variables to right 
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Table 1 
Cluster Sample 









band 




cluster 


rcdshift 


(lO 44 /^ 2 erg s- 1 ) 


(keV) 


references- z; 


MS 1137.5 + 6625 


0.784 


5.4 


0.3-3.5 


D99;GL94 


MS 0451.6- 0305 


0.550 


20.0 


0.3-3.5 


GL94;GL94 


CI 0016 + 16 


0.546 


14.6 


0.3-3.5 


DG92;GL94 


RX J1347.5- 1145 


0.451 


73.0 


0.1 - 2.4 


S95;S97 


Abcll 370 


0.374 


11. 7 b 


0.1 - 2.4 


M88;AE99 


MS 1358.4+ 6245 


0.327 


10.6 


0.3-3.5 


GL94;GL94 


Abell 1995 


0.322 


13.4 


0.1 - 2.4 


P00;B00 


Abell 611 


0.288 


8.6 


0.1 - 2.4 


C95;B00 


Abell 697 


0.282 


19.2 


0.1 - 2.4 


C95;B00,E98 


Abell 1835 


0.252 


32.6 


0.1 - 2.4 


SR99;B00,E98 


Abell 2261 


0.224 


20.6 


0.1 - 2.4 


C95;B00 


Abell 773 


0.216 


12.1 


0.1 - 2.4 


SR99;B00,E98 


Abell 2163 


0.202 


37.5 


0.1 - 2.4 


SR99;E96 


Abell 520 


0.202 


14.5 


0.1 - 2.4 


GL94;B00,E98 


Abell 1689 


0.183 


20.7 


0.1 - 2.4 


SR91;E96 


Abell 665 


0.182 


15.7 


0.1 - 2.4 


SR91;B00 


Abell 2218 


0.171 


8.2 


0.1 - 2.4 


L92;B00 


Abell 1413 


0.142 


10.9 


0.1 - 2.4 


SR99;B00,E98 


a Computed for a 


flat CIm = 1 


universe. 







b Converted the 2-10 keV flux in AE99 to the 0.1-2.4 keV band (approximate 
factor of 0.9 determined from cooling function calculation). 



REF: AE99-Arnaud & Evrard 1999; B00-B6hringer et al. 2000; C95-Crawford et al. 
1995; D99-Donahue et al. 1999; DG92-Dressler & Gunn 1992; E96-Ebeling et al. 
1996b; E98-Ebeling et al. 1998; GL94-Gioia & Luppino 1994; L92-Le Borgne et al. 
1992; M88-Mellier et al. 1988; POO-Patel et al. 2000; S95-Schindler et al. 1995; 
S97-Schindler et al. 1997; SR91-Struble & Rood 1991; SR99-Struble & Rood 1999; 
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ascension and declination are commonly called u and v, respectively, and the Fourier domain is commonly referred to as 
the u-v plane. The real and imaginary Fourier component pairs as a function of u and v are called visibilities. 

The finite size of each telescope dish results in an almost Gaussian response pattern, known as the primary beam. The 
product of the primary beam and the sky brightness distribution is equivalent to a convolution in the Fourier domain. 
The primary beams are measured using holography data for both OVRO and BIMA. The main lobe of the primary beams 
are well fit by a Gaussian with a full width at half maximum (FWHM) of 4'.2 for OVRO and 6'.6 for BIMA at 28.5 GHz. 
However, we use the measured primary beam profiles for our analysis. 

The primary beam sets the field of view. The effective resolution, called the synthesized beam, depends on the sampling 
of the u-v plane and is therefore a function of the configuration of the telescopes and the declination of the source. The 
cluster SZE signal is largest on the shortest baselines (largest angular scales). The shortest possible baseline is set by the 
diameter of the telescopes, D. Thus the system is not sensitive to angular scales larger than about X/2D, which is ~ 2'.8 
for BIMA observations and ~ 1'.7 for OVRO observations. The compact configuration used for our observations yields 
significant SZE signal at these angular scales, but filters out signal on larger angular scales. Because of the spatial filtering 
by the interferometer, it is necessary to fit models directly to the data in the u-v plane, rather than to the deconvolved 
image. 

Interferometers simultaneously measure both the cluster signal and the point sources in the field. The SZE signal is 
primarily present in the short baseline data while the response of an interferometer to a point source is independent of 
the baseline. Therefore, observations with a range of baselines allow us to separate the extended cluster emission from 
point source emission. We show an example of this after first presenting details of deconvolved 30 GHz images (see Jones 
et al. 1993; Carlstrom et al. 2000, for additional examples). 

3.1.2. Data Reduction 

The data are reduced using the MIRIAD (Sault et al. 1995) software package at BIMA and using MMA (Scovillc 
et al. 1993) at OVRO. In both cases, data are excised when one telescope is shadowed by another, when cluster data are 
not straddled by two phase calibrators, when there are anomalous changes in instrumental response between calibrator 
observations, or when there is spurious correlation. For absolute flux calibration, we use observations of Mars and adopt 
the brightness temperature from the Rudy (1987) Mars model. For observations not containing Mars, calibrators in those 
fields are bootstrapped back to the nearest Mars calibration (see Grego et al. 2001 for more details). The observations 
of the phase calibrators over each summer give us a summer-long calibration of the gains of the BIMA and OVRO 
interferometers. They both show very little gain variation, changing by less than 1% over a many-hour track, and the 
average gains remain stable from day to day. In fact, the gains are stable at the ~ 1% level over a period of months. 

Table 2 
Cluster Data 



Interferometric SZE Data ROSAT Data 







OVRO (hr) 






BIMA (hr) 




PSPC 


HRI 


cluster 


1994 


1995 


1996 


1998 


1996 


1997 


1998 


2000 


(ks) 


(ks) 


MS1137 












40 


48 






99.1 


MS0451 






30 












15.4 


45.9 


C10016 


87 


13 






29 a 


8 






41.6 


70.2 


R1347 
















20.0 




36.1 


A370 






33 






26 








31.9 


MS1358 






9 


7 






70 




22.1 


29.2 


A1995 






58 








50 






37.6 


A611 




45 


12 














17.2 


A697 




47 
















27.8 


A1835 














27 




8.5 


2.8 


A2261 




40 










3 






16.1 


A773 


57 


9 










18 






16.5 


A2163 






25 


12 


2 


10 


11 




11.7 


35.8 


A520 








7 


13 a 


23 


20 




4.7 


12.6 


A1689 




26 








16 






13.5 


22.5 


A665 










38 a 


24 






37.0 


98.3 


A2218 




64 


6 




20 a 


12 






42.5 


35.5 


A1413 










11 


17 






7.5 


18.6 



a Contains 1996 BIMA data with delay loss problem; data only used to make images and 
not in the analysis. 
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3.1.3. Data Visualization: 30 GHz Images 

Here we present deconvolved images of our 30 GHz intcrferomctric observations. However, we stress that these images are 
made to demonstrate the data quality. The model fitting is performed in the Fourier plane, where the noise characteristics 
of the data and the spatial filtering of the interferometer are well understood. The SZE and X-ray image overlays of 
Figure 2 show that the region of the cluster sampled by the interferometric SZE observations and the X-ray observations 
is similar for the clusters in our sample. In addition, the interferometer measures a range of angular scales, which is 
not apparent from the images in Figure 2. Images showing examples of our SZE data at varying resolutions appear in 
Carlstrom et al. (1996) for C10016 and Carlstrom et al. (2000) for R1347. 

Point sources are identified from SZE images created with DIFMAP (Pearson ct al. 1994) using only the long baseline 
data (> 2000 A) and natural weighting (cr -2 weight). Approximate positions and fluxes for each point source are obtained 
from this image and used as inputs for the model fitting discussed in §4.2. The data are separated by observatory, 
frequency, and year to allow for temporal and spectral variability of the point source flux. The positions and fluxes of 
the detected point sources from the model fitting are summarized in Table 3. Also listed are the corresponding 1.4 GHz 
fluxes for these sources from the NRAO VLA Sky Survey (NVSS) (Condon et al. 1998) and the 5 GHz and 15 GHz fluxes 
for sources in the three cluster fields surveyed by Moffet & Birkinshaw (1989). The uncertainty in the positions of the 
point sources is ~ ±3" at 68.3% confidence based on model fits of point sources described in §4.2. Figure 1 shows the 
30 GHz high resolution (> 2000 A) maps (color scale) with NVSS 1.4 GHz contours. The color scale wedge above each 
image shows the range in the map in units of mJy beam -1 . Contours arc multiples of twice the NVSS rms (rms ~ 0.45 
mJy beam -1 ). The FWHM of the 30 GHz synthesized beam is shown in the lower left hand corner of each panel and the 
45" FWHM beam of the NVSS survey is shown in the lower right hand corner of each panel. Table 4 summarizes the 
sensitivity and the FWHM of the synthesized beam of the high resolution maps used to find point sources in each field 
and shown in Figure 1. 

Figure 2 shows the deconvolved SZE image contours overlaid on the X-ray images for each cluster in our sample. 
Negative contours are shown as solid lines and the contours are multiples of twice the rms of each image. The images 
for MS0451 and C10016 have been published previously (Reese et al. 2000), but we include them here so that the entire 
sample appears together. We use DIFMAP (Pearson et al. 1994) to produce the naturally weighted SZE images. If any 
point sources are detected in the cluster field, they are subtracted from the data and a Gaussian taper applied to the 
visibilities to emphasize brightness variations on cluster scales before the image is deconvolved (CLEANed). The half- 
power radius of the Gaussian taper applied varies between 1000 A and 2000 A, depending of the observatory and telescope 
configurations used during the observations. Typically a 1000 A half-power radius taper is applied to BIMA data and a 
2000 A half-power radius taper is applied to OVRO data. The full-width at half maximum (FWHM) of the synthesized 
(restoring) beam is shown in the lower left hand corner of each image. Table 4 summarizes the rms sensitivities and the 
FWHM's of the synthesized beams of the tapered maps shown in Figure 2, as well as the corresponding statistics for the 
high resolution (> 2000 A) images. In addition, Table 4 lists the Rayleigh- Jeans (RJ) brightness sensitivities for each 
tapered, deconvolved image. 

3.1.4. Point Source Identification Using Spatial Filtering 

The identification and removal of point sources by taking advantage of the spatial filtering of the interferometer is 
illustrated in the panels of Figure 3 for the BIMA SZE C10016 data. During the maximum likelihood joint analysis (see 
§4.2), radio point sources identified from this procedure are modeled and fit for directly in the u-v plane. Each panel covers 
the same angular region, roughly 20' on a side, and each panel shows the FWHM of the synthesized beam in the lower left 
hand corner. Above each image is the color scale mapping showing the flux density of the map in units of mJy beam -1 . 
Panel a shows the "natural" image, which includes all of the data. There is smooth, extended, negative emission in the 
center of the map; this is the SZE decrement of C10016. There is also a bright spot roughly 2' south of the cluster that 
may be a point source. The large scale symmetric pattern is the synthesized beam of the low resolution data (compare 
to panel c); even when all baselines are considered, the SZE signal dominates. Figure 36 shows the high resolution map 
using data with projected baselines > 2000 A only. The point source shows up easily now with the characteristic shape 
of the synthesized beam for these data. We remove the point source by CLEANing. A Gaussian u-v taper (half-power 
radius of 1000 A) is then applied to the full data set to emphasize the short baselines, corresponding to the angular scales 
typical of galaxy clusters, and shown in Figure 3c. The cluster is apparent as is the symmetric pattern of the synthesized 
beam. Deconvolving (CLEANing) the tapered image results in panel d, which appeared in Figure 2 overlaid on X-ray 
data. The contours are multiples of twice the rms of the map. 

3.2. ROSAT X-ray Data 

We use archival Rontgen Satellite (ROSAT) data from both the Position Sensitive Proportional Counter (PSPC) and 
High-Resolution Imager (HRI) instruments. The live times of the observations we use are listed in Table 2 for both PSPC 
and HRI observations. 

3.2.1. Data Reduction 

We use the Snowden Extended Source Analysis Software (ESAS) (Snowden et al. 1994; Snowden 1998) to reduce the 
data. We use the ESAS software to generate a raw counts image, a noncosmic background image, and an exposure map 
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Table 3 
Radio Point Sources 



R.A. a (J2000) 


Dec. a (J2000) 
(d ' ") ' 


^30.0 


F28.5 


Fi5 b 


F 5 h 


F 1A C 


Field (h m s) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 



MS1137 

MS0451 

C10016 

R1347 

A370 

MS1358 

A1995 

A611 

A697 

A1835 

A2261 
A773 



04 


54 


22 


1 


-03 


01 


25 


00 


18 


31 


1 


+ 16 


20 


45 


13 


47 


30 


7 


-11 


45 


09 


02 


39 


55 


5 


-01 


,34 


06 


13 


59 


50 


G 


+62 


31 


05 


14 


53 


00 


5 


+58 


03 


19 


14 


01 


02 





+02 


52 


42 


14 


01 


00 


5 


+02 


51 


53 


17 


22 


17 


1 


+32 


09 


14 



1.41 

0.84 
0.58 



-0.26 
-0.26 



+0.09 
0.09 

+0.05 
-0.05 



A2163 


16 


15 


43.3 


-06 


08 


40 


1 44+012 

1 -^-0.12 


A520 


04 


54 


01.1 


+02 


57 


47 






04 


54 


17.0 


+02 


55 


32 






04 


54 


20.3 


+02 


54 


56 




A1689 


13 


11 


31.6 


-01 


19 


33 


1 oo+O.lO 
1 - oo -0.10 




13 


11 


30.1 


-01 


20 


37 


u -*°-0.09 


A665 


08 


31 


30.9 


+65 


52 


35 




A2218 


16 


35 


22.1 


+66 


13 


23 


4 29+ ' 21 

^■ zy -0.21 




16 


35 


47.7 


+66 


14 


4G 


1 36 +0 10 




16 


3G 


16.0 


+66 


14 


23 


2 41 + ' 29 

z -* i -0.29 


A1413 


11 


55 


08.7 


+23 


2G 


17 





1.86 
9.11 
10.81 
0.77 
1.61 
0.58 



0.26 
0.26 
+ 1.97 
1.97 
+0.19 
0.19 
+0.07 
-0.07 
+0.17 
-0.17 
0.04 
0.04 



+0.14 
0.14 
+0.15 
J -0.15 

10.801°;* 



2.76 
1.16 



+0.08 
-0.08 
+0.23 
0.23 
+0.10 



1.44 
7.97 

l-01-o.io 
1 03+ 012 

1 - UO -0.12 

1 r-,+0.09 

1 - oi -0.09 

0.42lg;g 
4.83±S;H 

4 4Q+0-20 
^■^°-0.20 

1 59+ 011 

llOy -0.11 

,+0.30 
0.30 
+0.23 
-0.23 



3.13 
2.01 



25 



+1.5 
-1.5 



1 9 7+0.3 

lz -'-0.3 



-,+0.6 
-0.6 
(+0.6 
l -0.6 

2 +0 - 5 

z - u -0.5 



5.0] 
2.4" 1 



84 



+1.1 

-1.1 



25. 
2. 
3. 
4. 



7+1.6 
' -1.6 
0+0.2 

S -0.2 
7+0.3 
' -0.3 
9+0.1 
z -0.1 



14.918;? 
269.3l|i 

47.6li;S 
7+1.1 
-1.1 



11.7 



3+0.9 

="-0.9 



41.4 



+1.9 
-1.9 



24.31J; 



6.7 
15.3 



+0.5 
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a Positions from SZE observations. 
b From Moffet & Birkinshaw (1989). 
c From Condon et al. (1998). 
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Fig. l. — -High resolution (> 2000 A) 30 GHz (color scale) with NVSS 1.4 GHz contours. The color scale wedge above 
each image shows the range in the 30 GHz map in units of mJy beam -1 . Contours are multiples of twice the NVSS rms 
of ~ 0.45 mJy beam -1 . The FWHM of the 30 GHz synthesized beam is shown in the lower left hand corner of each panel 
and the 45" FWHM beam of the NVSS survey is shown in the lower right hand corner of each panel. The 30 GHz image 
statistics are summarized in Table 4. 
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Fig. 2. — -SZE (contours) and X-ray (color scale) images of each cluster in our sample. Negative contours are shown as 
solid lines. The contours are multiples of 2 a and the FWHM of the synthesized beams are shown in the bottom left 
corner. The X-ray color scale images are raw counts images smoothed with Gaussians with a = 15" for PSPC data and 
a = 5" for HRI data. There is a color scale mapping for the counts above each image. The 30 GHz image statistics are 
summarized in Table 4. 
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Fig. 3. — -Panels illustrating how the spatial filtering of the interferometer is used to disentangle the point source emission 
from the SZE emission. Each panel is roughly 20' on a side with the FWHM of the synthesized beam shown in the lower 
left hand corner. The mapping of the color scale is shown above each panel in units of mJy beam -1 . Panels: a) "natural" ; 
b) high resolution (> 2000 A); c) tapered map after point source removal; and d) deconvolved SZE image (contours are 
multiples of 2 a). See text for details. 
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for the HRI (0.1-2.4 keV) data and for each of the Snowden bands R4-R7 (PI channels 52 — 201; approximately 0.5 — 2.0 
keV) for the PSPC data, using a master veto rate (a measure of the cosmic-ray and 7-ray backgrounds) of 200 counts s _1 
for the PSPC data. We examine the light curve data of both instruments looking for time intervals with anomalously high 
count rates (short-term enhancements) and for periods of high scattered solar X-ray contamination. Contaminated and 
anomalously high count rate data are excised. The Snowden software produces 512 x 512 pixel images with 14'.'947 pixels 
for the PSPC and 5'.'0 pixels for the HRI. For the PSPC, final images for all of the R4-R7 bands together are generated by 
adding the raw counts images and the background images. Each Snowden band has a slightly different effective exposure 
map and there is an energy dependence in the point spread function (PSF). Thus, we generate a single exposure image 
and a single PSF image by combining cluster photon-weighted averages of the four exposure images and the four PROS 
(Worrall et al. 1992; Conroy et al. 1993) generated on-axis PSF images. The cluster photon-weighting is determined 
using the background subtracted detected photons within a circular region centered on the cluster. The region selected 
to construct the weights is the largest circular region encompassing the cluster which contains no bright point sources, 
typically a 15 pixel radius. 

3.2.2. X-ray Images and Data Properties 

We show smoothed X-ray raw counts images in Figure 2 (color) with SZE image contours overlaid. PSPC images 
are shown when available and HRI images otherwise. Table 2 summarizes the on-source integration time of the ROSAT 
observations of the clusters in our sample for both the PSPC and HRI. These images roughly contain a few thousand 
cluster counts. PSPC images are smoothed with Gaussians with a = 15" and HRI images with a = 5". The color scale 
wedge above each figure shows the mapping between color and detector counts. 

3.2.3. X-ray Spectral Data 

We used published temperatures, metallicities, and H I column densities from observations with the Advanced Satellite 
for Cosmology and Astrophysics (ASCA). Temperatures and metallicities for most of the clusters in our sample appear in 
Allen & Fabian (1998a), Allen & Fabian (1998b), and Allen (2000). When there is a detailed account of the analysis for a 
particular cluster we use those results instead. When fitted metallicities are unavailable, we adopt a 0.2 solar metallicity 
with a 100% uncertainty. We use fitted H I column densities when available, otherwise those from 21 cm surveys of our 
galaxy (Dickey & Lockman 1990) are adopted. We assign a conservative 50% uncertainty to the column densities adopted 



Table 4 
SZE Image Statistics 



Tapered a High Resolution' 3 







a 


Beam 






a 


Beam 


Cluster 


Observatory 


(/xjy beam -1 ) 


(arcsec) 


( M K) 


Observatory 


(/^Jy beam -1 ) 


(arcsec) 


MS1137 


BIMA 


120 


90 


X 


94 


21 


BIMA 


105 


20 x 24 


MS0451 


OVRO 


90 


44 


X 


69 


45 


OVRO 


60 


19 x 25 


C10016 


BIMA 


250 


81 


X 


101 


46 


BIMA 


220 


14 x 30 


R1347 


BIMA 


307 c 


93 


X 


94 c 


53 c 


BIMA 


245 


17 x 27 


A370 


OVRO 


60 


5G 


X 


80 


19 


OVRO 


70 


17 x 23 


MS1358 


BIMA 


140 


9G 


X 


98 


22 


BIMA 


120 


17 x 20 


A1995 


BIMA 


134 


70 


X 


77 


37 


OVRO 


65 


17 x 20 


A611 


OVRO 


60 


48 


X 


58 


32 


OVRO 


45 


20 x 39 


A697 


OVRO 


65 


50 


X 


53 


37 


OVRO 


50 


19 x 33 


A1835 


BIMA 


213 


87 


X 


121 


30 


BIMA 


190 


18 x 22 


A2261 


OVRO 


85 


49 


X 


53 


49 


OVRO 


75 


19 x 35 


A773 


BIMA 


260 


91 


X 


99 


43 


OVRO 


90 


19 x 26 


A2163 


BIMA 


300 


90 


X 


104 


48 


OVRO 


85 


19 x 30 


A520 


BIMA 


180 


90 


X 


101 


30 


OVRO 


80 


12 x 16 


A1689 


BIMA 


320 


93 


X 


94 


55 


OVRO 


72 


18 x 49 


A665 


BIMA 


160 


93 


X 


99 


26 


BIMA 


150 


16 x 26 


A2218 


BIMA 


200 


93 


X 


99 


31 


OVRO 


50 


21 x 22 


A1413 


BIMA 


250 


93 


X 


99 


44 


BIMA 


210 


16 x 26 



a Gaussian taper with FWHM of 2000 A for OVRO data and 1000 A for BIMA data. 
b Using only data with ^J(u 2 + v 2 ) > 2000 A. 
c Used Gaussian taper with FWHM of 1500 A. 
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from 21 cm surveys of our galaxy. Table 5 summarizes our adopted electron temperatures, metallicities, and column 
densities with references for the sources of this information. 

Temperatures for many of our clusters also appear in Mushotzky & Scharf (1997). Multiple temperature determinations 
agree within the 1 a intervals for most of the clusters in our sample. The measurements overlap within 2 a in the worst 
cases, i.e., for the clusters, MS1358, A1995, A2163, A1689, and A1413. The temperatures, metallicity, column densities, 
and redshift of the cluster are used to determine the X-ray cooling functions and the conversion factor between detector 
counts and cgs units, E (see §3.2.4 below for details). The cooling functions and conversion factors are summarized in 
Table 6. 

3.2.4. X-ray Cooling Function 

The X-ray cooling function enters the distance calculation linearly and indirectly as the conversion between detector 
counts and cgs units (see §4.1). We use a Raymond-Smith (1977) spectrum to describe the hot ICM, which includes 
contributions from electron- ion thermal bremsstrahlung, line emission, recombination, and two photon processes. We 
replace the non-relativistic bremsstrahlung calculation in the Raymond-Smith model with the relativistic calculation of 
Gould (1980). A discussion of this calculation appears in Reese et al. (2000). 

The cooling function results for the ROSAT data used in our analysis arc summarized in Table 6, where A eH o is the 
cooling function in cgs units, Af^g is the cooling function in detector units, E is the conversion between counts and cgs 
units, and A bol is the bolometric cooling function. The cooling functions with relativistic corrections are typically 1.05 
times the Raymond-Smith "uncorrected" value for the clusters in our sample. 



Table 5 
X-ray Spectral Information 





kT e 








cluster 


(keV) 


[Fc/H] 


(xlO 20 cm- 2 ) 


ref-T e ; [Fe/H]; N H 


MS1137 


K 7+1.3 
°-'-0.7 


0.4318;^ 


1.00 a 


D99;D99;D99 


MS0451 


10 4+ 10 


1 5 +a07 


o oo+ - 40 


D96;D96;D96 


C10016 


7 cc+0.72 
'■ oo -0.58 


07+ 11 

u - u ' -0.07 


r. KO+0-41 
o - o9 -0.36 


HB98;HB98;HB98 


R1347 


O Q+0-7 
a -°-0.6 


n oq+0.06 
u -°°-0.06 


io.ooH-88 


S97;S97;S97 


A370 


°-°-0.5 


3+ 01 


3.1 b 


098;098;G 


MS1358 


7.48±° j° 


n OO+0.09 

u - oz -0.09 


1.93 b 


AF98;AF98b;G 


A1995 


8.5918;?? 


14+ 007 

u - 1 ^-0.07 


5-Oli-j 


P00;P00;P00 


A611 


6.612:8 


0.20 c 


4.99 b 


HPC;-;G 


A697 


9.81°;? 


0.20 c 


3.41 b 


HPC;-;G 


A1835 


o 01 +0.19 


S5+ 004 

u -°°-0.03 


2.32 b 


AF98;AF98b;G 


A2261 


o no+0.37 
8 - Sz -0.32 


n 09+O.O6 
u - oz -0.05 


3.28 b 


AF98;AF98b;G 


A773 


O oq+0.41 


21 + 05 


1.44 b 


AF98;AF98b;G 


A2163 


i o 2+ 11 

iz ' z -0.7 


40+ 09 

u -^ u -0.08 


16.501?;?° 


M96;EAB95;EAB95 


A520 


o oq+0.46 


n 14+0.06 


7.80 b 


AF98;AF98b;G 


A1689 


9 fifi +0 22 

y.oo„ 20 


u - zy -0.03 


1.82 b 


AF98;AF98b;G 


A665 


9 D3+ - 35 


n 99 +0.04 
u - zz -0.05 


4.24 b 


AF98,AF98b;G 


A2218 


7 05 +a22 

' - U0 -0.21 


n 1s +0.04 
U - iS -0.04 


3.24 b 


AF98;AF98b;G 


A1413 


7 c 4 +o.i7 


n 9 fi +0.03 
U - zS -0.03 


2.19 b 


AF98;AF98b;G 



a Adopted value (see ref D99) and adopted 50% uncertainty. 

b Galactic value (Dickey & Lockman 1990) with adopted 50% uncertainty. 

c Adopted value with assumed 100% uncertainty. 

NOTE: Uncertainties are 68% confidence. 
REF: AF98-Allen & Fabian 1998a; AF98b-Allen & Fabian 1998b; D96- 
Donahue 1996; D99-Donahue et al. 1999; EAB95-Elbaz et al. 1995; G- 
Dickcy & Lockman 1990; HPC-Hughes, private communication; HB98- 
Hughes & Birkinshaw 1998; M96-Markevitch et al. 1996; 098-Ota et al. 
1998; POO-Patcl et al. 2000; S97-Schindlcr et al. 1997 
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Table 6 
X-ray Cooling Functions 



PSPC HRI 



cluster 


A eH o a 


A det b 
yi eH0 






A det b 
Yi eH0 




S c 


n e 
n H 


Mh 

Me 


A d 


MS1137 


7.751 










1.765 


2 


,461 


1 


.202 


2.146 


MS0451 


6.948 


3 


263 


1 


.373 


1.470 


3 


,050 


1 


.198 


2.702 


C10016 


6.922 


3 


003 


1 


.489 


1.289 


3 


,471 


1 


.197 


2.260 


R1347 


6.922 










1.167 


4 


,089 


1 


.201 


2.643 


A370 


6.790 










1.627 


3 


,037 


1 


.200 


2.223 


MS1358 


6.717 


3 


,785 


1 


.336 


1.793 


2 


,821 


1 


.200 


2.371 


A1995 


6.434 










1.434 


3 


395 


1 


.198 


2.448 


A611 


6.511 










1.489 


3 


395 


1 


.199 


2.175 


A697 


6.334 










1.570 


3 


,148 


1 


.199 


2.647 


A1835 


6.462 


3 


.839 


1 


.344 


1.775 


2 


,909 


1 


.201 


2.496 


A2261 


6.359 










1.649 


3 


,150 


1 


.200 


2.570 


A773 


6.171 










1.870 


2 


,715 


1 


.199 


2.581 


A2163 


6.135 


2 


.555 


1 


.998 


1.021 


5 


,000 


1 


.202 


3.064 


A520 


6.119 


3 


.209 


1 


.586 


1.337 


3 


,807 


1 


.198 


2.411 


A1689 


6.158 


3 


,899 


1 


.335 


1.835 


2 


,836 


1 


.200 


2.673 


A665 


6.102 


3 


616 


1 


.428 


1.570 


3 


,289 


1 


.199 


2.550 


A2218 


6.112 


3 


,789 


1 


.378 


1.691 


3 


,086 


1 


.198 


2.237 


A1413 


6.133 


4 


000 


1 


.343 


1.856 


2 


,894 


1 


.200 


2.361 



a Units are xlO -24 erg s _1 cm 3 . The emissivity in the cluster frame 
integrated over the ROSAT band (0.5-2.0 keV) redshiftcd to the cluster 
frame. 

b Units are xlO -13 cnts s _1 cm 5 . The emissivity in the detector frame 
accounting for the response of the instrument. 

c Units are x 

10 _n erg s -' cm- _ Th{j 

conversion of detector units to 

cnt s 

cgs units including the (1 + z) factor between energy and counts. 
d Units are xlO -23 erg s _1 cm 3 . The bolometric emissivity. 
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4. METHOD 

4.1. Angular Diameter Distance Calculation 

The calculation begins by constructing a model for the cluster gas distribution. We use a spherical isothermal f3 model 
to describe the ICM. With this model, the cluster's extent along the line of sight is the same as that in the plane of 
the sky. This is clearly invalid in the presence of cluster asphericities. Thus cluster geometry introduces an important 
uncertainty in SZE and X-ray derived distances. In general, clusters are dynamically young, are aspherical, and rarely 
exhibit projected gas distributions which are circular on the sky (Mohr ct al. 1995). We currently cannot disentangle 
the complicated cluster structure and projection effects, but numerical simulations provide a good base for understanding 
these difficulties. The effects of asphcricity contribute significantly to the distance uncertainty for each cluster, but are 
not believed to result in any significant bias in the Hubble parameter derived from a large sample of clusters (Sulkancn 
1999). 

The spherical isothermal (3 model has the form (Cavaliere & Fusco-Femiano 1976, 1978) 

9 v -3/3/2 



n e (r) = n e0 I 1 + 1 , (1) 

where n e is the electron number density, r is the radius from the center of the cluster, r c is the core radius of the ICM, 
and (3 is the power law index. With this model, the SZE signal is 

C k T ( 9 2 \ (1 " 3/3) / 2 

AT = f {xTc) T CMB D A Jd(a T n e -^ = AT Q (l + ^ J , (2) 

where AT is the thermodynamic SZE temperature decrement /increment, /(x,r e ) is the frequency dependence of the 
SZE with x = his/kT CMB , T CMB (=2.728 K; Fixsen et al. 1996) is the temperature of the CMB radiation, k B is the 
Boltzmann constant, cr T is the Thompson cross section, m e is the mass of the electron, c is the speed of light, ATq is the 
central thermodynamic SZE temperature decrement/increment, 9 is the angular radius in the plane of the sky and 9 C the 
corresponding angular core radius, and the integration is along the line of sight I = D A (. The frequency dependence of 
the thermal SZE is 

/(x,T e) = (zj^j - 4^ (1 + S SZE (x,T e )) (3) 

where 5 SZE (x,T e ) is the relativistic correction to the frequency dependence. In the non-rclativistic and Rayleigh- Jeans 
(RJ) limits, f( x ,T e ) ~ * — 2. We apply the relativistic corrections 8 SZE (x,T e ) to fifth order in kT e /m e c 2 (Itoh et al. 1998), 
which agrees with other works (Rephaeli 1995; Rephaeli & Yankovitch 1997; Stebbins 1997; Challinor & Lasenby 1998; 
Sazonov & Sunyaev 1998a, b; Molnar & Birkinshaw 1999; Dolgov et al. 2001) for clusters with k B T e < 15 keV, satisfied by 
all the clusters in our sample. This correction decreases the magnitude of fr x ,T e ) by < 5% (typically 3%) for the clusters 
considered here. 

The X-ray surface brightness is 

If / 02\(l-60)/2 

S - = 4 n (l + z )4 DA J dCnen » Ae » =S ^[ 1+ 0l) ' W 

where S x is the X-ray surface brightness in cgs units (erg s _1 cm~ 2 arcmin~ 2 ), z is the redshift of the cluster, n H is the 
hydrogen number density of the ICM, A eH = A eH (T e , abundance) is the X-ray cooling function of the ICM in the cluster 
rest frame in cgs units (erg cm 3 s _1 ) integrated over the redshifted ROSAT band, and S x q is the X-ray surface brightness 
in cgs units at the center of the cluster. Since the X-ray observations are in instrument counts, we also need the conversion 
factor between detector counts and cgs units, S (S x o = S x qH), discussed in detail in Reese et al. (2000) along with a 
description of the calculation of A eH , which includes relativistic corrections (Gould 1980) to the Raymond & Smith (1977) 
spectrum. The normalizations, ATq and S x o, used in the fit include all of the physical parameters and geometric terms 
that come from the integration of the (3 model along the line of sight. 

One can solve for the angular diameter distance by eliminating n e o (noting that n H = n e fi e / /x H where nj = p/ /J,jm p for 
species j) yielding 



n (AT ) 2 /m e c 2 V A 



, 2 



r(3/3 - 1) 
r(3/3) 



S x0 \k B T e0 J 47r 3 / 2 / ( 2 X) T e) 7 1 c 2 MB a2(l + z) 

where T(x) is the Gamma function. Similarly, one can eliminate D A instead and solve for the central density n e Q. 
More generally, the angular diameter distance is 

/ \ 2 A 

(AT ) 2 f m e c 2 \ 2 A eH0 [i e /p H 1 / (n^) 7^ dr] l*=o 

t k B T e0 J Air ff xTe) T 2 MB a 2 (1 + zf 9, 



(5) 



„ V"^u; '"e- \ "chupc; m ^ - y-eu/ £^eH0 " — 

- J-2 T2 ^2 f1 x ,\4 fl r n y 



where 9 C is the characteristic angular scale of the galaxy cluster whose exact meaning depends on the ICM model (the 
core radius for the isothermal (3 model) and 77 = Q/9 C = £/r c is the line of sight length in units of the characteristic radius, 
r c = 9 C D A . For simplicity in notation, we have assumed that the density and temperature models are normalized at the 
central value (denoted with 0) though any location for the normalization is allowed. The above integrals are along the 
central line of sight, denoted as zero projected radius R = 0. The T functions and the factor of 7r i/2 in equation (5) come 
from the integration of the (3 model along the central line of sight for both the SZE and X-ray models. 
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4.2. Joint SZE and X-ray Model Fitting 

The SZE and X-ray emission both depend on the properties of the ICM, so a joint fit to the interferomctric SZE 
data and the PSPC and HR1 X-ray data provides the best constraints on those properties. Each data set is assigned a 
collection of parameterized models. Typically, SZE data sets are assigned a (3 model and point sources and X-ray images 
are assigned a (3 model and a cosmic X-ray background model. This set of models is combined for each data set to create 
a composite model which is then compared to the data. 

Our analysis procedure is described in detail in Reese et al. (2000). The philosophy behind the analysis is to keep 
the data in a reduced but "raw" state and run the model through the observing strategy to compare directly with the 
data. In particular, the interferometric SZE observations provide constraints in the Fourier (u-v) plane, so we perform 
our model fitting in the u-v plane, where the noise properties of the data and the spatial filtering of the interferometer 
are well defined. The SZE model is generated in the image plane, multiplied by the primary beam, and fast Fourier 
transformed to produce model visibilities. Wc then interpolate the model visibilities to the u and v of each data visibility 
and compute the Gaussian likelihood. For X-ray data, the model is convolved with the appropriate point spread function 
and the Poisson likelihood is computed pixel by pixel, ignoring the masked point source regions. 

Each data set is independent, and likelihoods from each data set can simply be multiplied together to construct the 
joint likelihood. Likelihood ratio tests can then be performed to get confidence regions or compare two models. Rather 
than working directly with likelihoods, C, wc work with S = — 21n(£). We then construct a Ax 2 -like statistic from 
the log likelihoods, AS = S n — S m i n where S ml „ is the minimum of the S function and S n is the S statistic where n 
parameters differ from the parameters at S m i n - The statistic AS is sometimes referred to as the Cash (1979) statistic 
and tends to a \ 2 distribution with n degrees of freedom for large n (Kendall & Stuart 1979, for example). This AS 
statistic is equivalent to the likelihood ratio test and is used to generate confidence regions and confidence intervals. For 
one interesting parameter, the 68.3% (~ la) confidence level corresponds to AS = 1.0. 

4.2.1. Model Fitting Uncertainty Estimation 

Uncertainties in the angular diameter distance from the fit parameters are calculated by varying the interesting pa- 
rameters to explore the AS* likelihood space. The most important parameters in this calculation are ATo, S x o, (3, and 
6 C . Radio point sources and the cosmic X-ray background affect ATq and S x o, respectively. As a compromise between 
precision and computation time, we systematically vary ATq, S x q, (3, and 9 C allowing the X-ray backgrounds for the PSPC 
and HRI to float independently while fixing the positions of the cluster (both SZE and X-ray), and the positions and flux 
densities of any radio point sources in the SZE cluster fields. We describe our estimation of the effects of point sources 
below. 

From this four dimensional AS hypcr-surface, we construct confidence intervals for each parameter individually as 
well as confidence intervals for D A due to S x o, AT , (i, and 9 C jointly. To compute the 68.3% confidence region we find 
the minimum and maximum values of the parameter within a AS of 1.0. We emphasize that these uncertainties are 
meaningful only within the context of the spherical isothermal j3 model. 

Measured Radio Point Sources 

Two methods of estimating the effect of the measured radio point sources in the cluster field arc examined, one which is 
reasonably quick and one that is more rigorous. For the quick estimate, we first determine the 1 a confidence limits on the 
flux density of each point source by varying the point source flux density while keeping the ICM parameters fixed at their 
best fit values. These are the uncertainties listed in Table 3, after correcting for the primary beam attenuation appropriate 
for each point source's distance from the pointing center. We then determine the change in the central decrement over the 
68.3% confidence region for the point source flux densities by fixing the point source flux density at the ±1 a values and 
varying ATo while fixing the ICM shape parameters at their best fit values. This is done for each point source in the field 
and all combinations of the ±1 a flux densities for fields with multiple point sources. We adopt the maximum percentage 
change in ATo as our uncertainty from radio point sources on the central decrement. The above procedure will be referred 
to as the quick estimate of the effects of measured radio point sources. We tested this estimate against marginalizing over 
the point source flux density by varying 6 C , (3, S x q, ATq, and point source flux for each point source, while fixing the X-ray 
background (simply saves computation time by isolating the point source flux issue). From the marginalized likelihood 
function we find the best fit and 68.3% uncertainty on AT and D A . The uncertainty from measured radio point sources, 
(T p t, is computed assuming the uncertainties add in quadrature from 

°pt = a mar ~ ° fix J CO 

where a mar and Ofi X are the uncertainties from the marginalized grids and the initial, point source fixed grids, respectively. 
Marginalizing over point source flux density was performed on two clusters with one point source each, A2261 and MS1358, 
and one cluster with two point sources, A1835. The quick estimation of the effects of point sources agrees to within 2% 
on D A (1% on ATo) with the marginalized likelihood analysis, just slightly over estimating the uncertainty due to point 
sources. The marginalization procedure is computationally intensive. Therefore, to save computation time, we use the 
quick procedure to estimate the effects of detected point sources on the central decrement. 

As an additional test, we explore the maximum likelihood parameter space (varying 9 C , (3, S x q, and ATo) with the 
point source fluxes fixed at the ±1 a values for our three test case clusters: A2261, MS1358, and A1835. The effects of 
point sources on the central decrement from this study agree within a few percent with both the quick and marginalized 
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procedures. This is what was originally done for MS0451 and C10016 (Reese et al. 2000), which is now shown to give 
essentially the same result as marginalizing over the point source flux. For all clusters, we use the updated, quick estimates 
of the effects of measured point sources. 

4.3. Model Fitting Results 

We apply the analysis procedure described above to all 18 of our galaxy clusters. The results from our maximum 
likelihood joint fit to the SZE and X-ray data are summarized in Table 7, which shows the best-fit ICM shape parameters 
and the uncertainties on each parameter from the model fit. 

So far, we have only shown the SZE data in the form of images though the data are recorded as visibilities. Figure 4 
shows the SZE u-v radial profiles for C10016 with a series of 3 panels illustrating the features of such profiles. These 
profiles are azimuthal averages in the Fourier plane plotted as a function of the radius in the u-v plane, \J u 2 + v 2 . The 
data are the points with error bars and the best fit (3 model from the joint SZE and X-ray analysis is the solid line 
averaged the same way as the data. The point sources are subtracted directly from the visibilities before constructing the 
u-v radial profiles. All of these panels are shown on the same scale for easy comparison. Also plotted are the residuals in 
units of the standard deviation, AV/<r=(data — modcl)/cr. For a circular cluster at the phase center (coincident with the 
pointing center), one expects a monotonic real component and a zero imaginary component. Clusters are rarely exactly 
centered at the phase/pointing center of our observations. Therefore, we shift the phase center to the cluster center before 
constructing the u-v radial profiles. The phase shifted radial profiles are shown in the upper panel of Figure 4 for both 
the real (left) and imaginary (right) components of the complex visibilities. The model provides a good fit to the data for 
a wide range of spatial frequencies. The middle panel shows the u-v radial profile when the phase center is not shifted 
to the center of the cluster. The off-center cluster introduces corrugation in the Fourier plane modifying the expected 
real component and introducing a non-zero imaginary component. In addition, asymmetry in the cluster will manifest 
itself as a non-zero imaginary component. Our model is symmetric and its imaginary component should be identically 
zero. The attenuation from the primary beam introduces asymmetry, producing a non-zero imaginary component. This is 
illustrated in the lower panel of Figure 4 showing the u-v radial profile including the phase center shift but not including 
the primary beam correction when computing the model. Notice the model is identically zero unlike the upper panel, 
where the asymmetry produced by the primary beam on the off-center cluster shows a small imaginary component. 

The u-v radial profiles for the real and imaginary components of the complex visibilities for each cluster in our sample 
are shown in Figure 5. Any point sources in the field are subtracted directly from the visibilities and the phase center 
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Fig. 4. — -The real and imaginary components of the complex visibilities plotted as a function of radius in the u-v plane 
for C10016. The points with error bars are the data and the best fit model from the joint SZE and X-ray analysis is the 
solid line, averaged the same way as the data. Residuals are shown in units of the standard deviation. Shown are the 
phase shifted and primary beam corrected (upper panel) , not phase shifted (middle panel) , and not corrected for primary 
beam (lower panel) versions. Not shifting the phase center to the center of the cluster (middle panel) shows an imaginary 
component from this offset. Not applying the primary beam attenuation to the model after shifting the phase to the 
center of the cluster (lower panel) shows the expected zero imaginary component; a real and symmetric image should 
have a real only Fourier transform. The asymmetry induced by the primary beam correction for the off-center cluster 
introduces a small imaginary component (sec upper panel). 
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Fig. 5. — -The real and imaginary components of the complex visibilities plotted as a function of radius in the u-u-planc. 
The points and error bars are the data and the solid line is the best fit model. Residuals are shown in units of the 
standard deviation. Any point sources in the cluster field have been subtracted directly from the visibilities and the phase 
center of the map has been shifted to the center of the cluster before making these radial averages. Non-zero signal in 
the imaginary components is due to the asymmetry in the cluster and the possible asymmetry introduced by the primary 
beam correction. The models provide good fits to the data for all the clusters in our sample. 
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is shifted to the center of the cluster before azimuthally averaging the real and imaginary components of the complex 
visibilities. The points with error bars are the data and the best fit model from the joint SZE and X-ray analysis is shown 
as a solid line, averaged the same way as the data. Also shown are the residuals in units of the standard deviation. A 
simple x 2 analysis of the SZE u-v radial profiles reveals that the models provide a good fit to the data for every cluster. 
The real and imaginary components are shown on the same scale for each cluster for easy comparison, though the scale 
on the residuals may change. The cluster with the most apparent imaginary component, A520, is also the cluster with its 
best fit center the furthest away from the pointing center, ~ 65". The primary beam attenuation introduces asymmetry 
and produces a non-zero imaginary component paralleled in the best fit model for A520. 

Figure 6 shows the X-ray radial surface brightness profiles and the best fit composite models for each cluster in the 
sample. Residuals in units of the standard deviation are also plotted, AS x /<7s T - A simple x 2 analysis of the radial profiles 
shows that, in general, the models provide a reasonable fit to the data over a large range of angular radii. There are 3 
clusters with residuals > 5 a in a few of the radial bins; A1835, A665, and A2281. In A1835 the model systematically 
undcrprcdicts the surface brightness for a few intermediate radial bins. This cluster contains a strong cooling flow (e.g., 
Peres et al. 1998), though shows no sign of a central emission excess over the best fit model. The other cooling flow clusters 
in our sample (see §6.1) do not exhibit large residuals in the X-ray model fit. Both A665 (Gomez et al. 2000; Kalloglyan 
et al. 1990; Geller & Beers 1982) and A2218 (Cannon et al. 1999; Girardi et al. 1997; Markevitch 1997; Kneib et al. 1995) 
exhibit complicated structure, possibly indicating a recent merger. Though the worst cases, the best fit models for these 
3 clusters still provide reasonable descriptions of the data. The only cluster that shows a marginally significant central 
X-ray surface brightness excess is A2218, which has not been identified with a cooling flow. We also note that the cooling 
flow clusters in our sample do not exhibit the largest residuals. 

5. DISTANCES AND THE HUBBLE CONSTANT 

We use the results from the maximum likelihood model fitting described in §4.2 and §4.3 to compute the angular 
diameter distance to each of our 18 galaxy clusters. Table 7 shows the derived angular diameter distances for each galaxy 
cluster as well as the best fit ICM shape parameters. The uncertainties on D A include the entire observational uncertainty 
budget, which arc shown for each cluster in Table 8. The uncertainties in the fitted parameters come from the procedure 
described in §4.2.1. 

The only other parameter that enters directly into the D A calculation is T e Q. Since D A tx T^q 2 , the uncertainty in D A 
due to T e o is listed as twice the fractional uncertainty on T e Q. The other parameters, column density and metallicity, as 
well as T e 0i affect the X-ray cooling function. We estimate the uncertainties in D A due to these parameters by taking 
their 68.3% ranges and seeing how much they affect the cooling function. The uncertainty on D A due to observations is 
dominated by the uncertainty in the electron temperature and the SZE central decrement. Note that changes of factors 
of two in metallicity result in a ~ 1% effect on D A . 

The column densities measured from the X-ray spectra are different from those from H I surveys (Dickey & Lockman 
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Fig. 6. — -X-ray surface brightness radial profiles with the best fit model and residuals in units of the standard deviation. 
In general, the models provide good fits to the data over a large range of radii. 
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+21.8 


-14.4 


-18.0 


-1.1 


- 2.2 


-1.0 


- 8.6 


-24.8 


+32.4 


+ 9.6 


+ 1.0 


+ 12.1 


+0.3 


+ 10.2 


+37.3 


-27.2 


-11.0 


-1.0 


-10.2 


-0.4 


-10.2 


-32.7 


+ 14.4 


+ 4.1 


+0.5 


+ 2.9 


+0.2 


+ 19.8 


+25.0 


-12.2 


- 4.6 


-0.5 


- 2.8 


-0.2 


-19.6 


-23.7 


+45.5 


+ 6.9 


+0.6 


+ 6.7 


+0.3 


+ 0.4 


+46.5 


-37.3 


- 7.8 


-0.8 


- 6.1 


-0.4 


- 0.6 


-38.5 


+30.3 


+ 6.0 


+0.7 


+ 5.2 


+0.2 


+ 10.0 


+32.9 


-31.7 


- 6.2 


-0.7 


- 4.8 


-0.2 


- 10.0 


-34.2 


+28.5 


+ 4.2 


+0.5 


+ 3.5 


+0.2 


+ 7.2 


+29.9 


-24.8 


- 4.5 


-0.5 


- 3.3 


-0.2 


- 7.2 


-26.4 



MS1137 

MS0451 

CL0016 

R1347 

A370 

MS1358 

A1995 

A611 

A697 

A1835 

A2261 

A773 

A2163 

A520 

A1689 

A665 

A2218 

A1413 



a Thc 68.3% uncertainties over the four-dimensional error surface for /3, 

J C1 

S x0 , and AT . 
h D A decreases as parameter increases. 
c Metallicity relative to solar. 
d Maximum effect from detected point sources. 
c Combincd in quadrature. 
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z (redshift) 

Fig. 7. — -SZE determined distances as a function of redshift. The error bars are 68.3% statistical uncertainties only. 
Also plotted are the theoretical angular diameter distance relations assuming Hq = 60 km s" 1 Mpc" 1 for three different 
cosmological models; the currently favored A cosmology fljif = 0.3, J7a = 0.7 (solid) cosmology; an open Ojvf = 0.3 
(dashed) universe; and a flat CIm = 1 (dotted) cosmology. 



1990). We use the column densities from X-ray spectral fits when possible since that includes contributions from non- 
neutral hydrogen and other elements which absorb X-rays. For MS0451 and C10016, using the survey derived column 
densities instead of the fitted values changes the angular diameter distance by ~ ±5% (Reese et al. 2000), which we 
include as a systematic uncertainty (see § 7). 

Figure 7 shows the SZE determined distances for each cluster as a function of redshift. Also plotted are the theoretical 
angular diameter distance relations assuming Hq = 60 km s _1 Mpc -1 for three different cosmologies: the currently favored 
A cosmology Hm = 0.3, Oa = 0.7 (solid) cosmology; an open Qm = 0.3 (dashed) universe; and a flat Qm = 1 (dotted) 
cosmology. The SZE distances arc beginning to probe the angular diameter distance relation. The uncertainties on D A 
in Figure 7 are the 68.3% statistical uncertainties only, including all of the statistical uncertainties in the calculation 
outlined above. We refer the reader to Carroll et al. (1992), Kolb & Turner (1990), and Peacock (1999) for derivations of 
the theoretical angular diameter distance relation. 

There is a known correlation between the (3 and 9 C parameters of the (3 model. Figure 8 illustrates this correlation and 
its effect on D A for MS1358, and A2261. The filled contours are the 1, 2, and 3 a AS confidence regions for f3 and 9 C 
jointly with the plus marking the best fit for each cluster. The lines are contours of constant D A in megaparsecs. With our 
intcrferometric SZE data, the contours of constant D A lie roughly parallel to the f3-6 c correlation, minimizing the effect 
of this correlation on the uncertainties of D A . Similar figures for MS0451 and C10016 appear in Reese et al. (2000), which 
show similar behavior. The alignment of the D A contours with the (3-6 c correlation is a general feature of our observing 
strategy. Different observing techniques will result in different behavior. Contours of constant D A have been found to be 
roughly orthogonal to this (3-9 c correlation for some single dish SZE observations (Birkinshaw & Hughes 1994; Birkinshaw 
et al. 1991). 

To determine the Hubble Constant, we perform a \ 2 fit to our calculated D A 's versus z for three different cosmologies. 
To estimate statistical uncertainties, we combine the uncertainties on D A listed in Table 8 in quadrature, which is only 
strictly valid for Gaussian distributions. This combined statistical uncertainty is symmetrized (averaged) and used in the 
fit. We find 

r 60+j km s- 1 Mpc" 1 ; ft M =0.3, ft A =0.7, 
H = < 56+4 km s- 1 Mpc" 1 ; O M =0.3, fi A =0.0, (8) 
{ 54+j km s- 1 Mpc" 1 ; f2 M =1.0, ft A =0.0, 

where the uncertainties are statistical only at 68.3% confidence. The statistical error comes from the \ 2 analysis and 
includes uncertainties from T e , the parameter fitting, metallicity, N H , and detected radio point sources (see Table 8). We 
have chosen three cosmologies encompassing the currently favored models. With this sample of clusters, there is a ~ 10% 
range in our inferred Hq due to the geometry of the universe. For the A cosmology, \ 2 = 16.5 with a corresponding 
reduced chi-squared of xl e d = 0-97. The difference in x 2 between the Qm = 0.3, il\ = 0.7 and the flat Qm = 1 universes 
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Fig. 8. — -Confidence regions from the joint SZE and X-ray fit for MS1358 and A2261. The filled regions are 1, 2, and 
3 a confidence regions for /3 and 9 C jointly (AS = 2.3, 6.2, 11.8), and the cross marks the best-fit (3 and 9 C . Solid lines 
are contours of angular diameter distance in megaparsecs. The D A contours lie roughly parallel to the j3-9 c correlation, 
minimizing the effect of this correlation on the uncertainties of D A . 

is roughly A\ 2 ~ 0.3, with the A cosmology having the lowest \ 2 - Clearly a larger sample of high redshift (z ~ 1) clusters 
is required for a determination of the geometry of the universe from SZE and X-ray determined direct distances to galaxy 
clusters (see §7). 

6. SOURCES OF POSSIBLE SYSTEMATIC UNCERTAINTY 

The absolute calibration of both the SZE observations and the PSPC and HRI directly affects the distance deter- 
minations. The absolute calibration of the interferometric observations is conservatively known to about 4% at 68.3% 
confidence, corresponding to a 8% uncertainty in H a (oc AT ~ 2 ). The effective areas of the PSPC and HRI are thought 
to be known to about 10%, introducing a 10% uncertainty into the Hq determination through the calculation of S. In 
addition to the absolute calibration uncertainty from the observations, there are possible sources of systematic uncertainty 

Table 9 

Hq Systematic Uncertainty Budget (%) 



Systematic Effect 



SZE calibration ± 8 

X-ray calibration ±10 

N H ±5 

Asphericity a ± 5 

Isothcrmality ±10 

Clumping —20 
Undetected radio sources 13 ±12 

Kinetic SZE a ± 2 

Primary CMB a < ±1 

Radio Halos — 3 

Primary Beam ± 3 

Total c til 



a Includes a 1/^/18 factor for our 
18 cluster sample. 

b Average of effect from the 18 clus- 
ter fields. 

c Combined in quadrature. 
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that depend on the physical state of the ICM and other sources that can contaminate the cluster SZE emission. Table 9 
summarizes the systematic uncertainties in the Hubble constant determined from our 18 cluster sample. 

6.1. Cluster Atmospheres and Morphology 
6.1.1. Asphericity 

Most clusters do not appear circular in radio, X-rays, or optical. Fitting a projected elliptical isothermal j3 model gives 
typical axial ratios that are close to the local average of 0.80 (Mohr et al. 1995). Under the assumption of axisymmetric 
clusters, the combined effect of cluster asphericity and its orientation on the sky conspires to introduce a ~ ±20% random 
uncertainty in Ho determined from one galaxy cluster (Hughes & Birkinshaw 1998). When one considers a large, unbiased 
sample of clusters, with random orientations, the errors due to imposing a spherical model are expected to cancel, resulting 
in a precise determination of H n . Recently, Sulkanen (1999) studied projection effects using triaxial (3 models. Fitting 
these with spherical models he found that the Hubble constant estimated from the sample was within 5% of the input 
value. We are in the process of using N-body and smoothed particle hydrodynamics (SPH) simulations of 48 clusters to 
quantify the effects of complex cluster structure on our results. 

A 20% effect from one cluster implies a 5% (~ 20/\/l8) effect for a sample of 18 clusters. Therefore, we include at 5% 
effect from asphericity for our cluster sample. 

6.1.2. Temperature Gradients 

Departures from isothermality in the cluster atmosphere may result in a large error in the distance determination from 
an isothermal analysis; moreover, an isothermal analysis of a large cluster sample could lead to systematic errors in the 
derived Hubble parameter if most clusters have similar departures from isothermality (Birkinshaw & Hughes 1994; Inagaki 
et al. 1995; Holzapfel et al. 1997b). The ROSAT band is fairly insensitive to temperature variations, showing a ~ 10% 
change in the PSPC count rate for a factor of 2 change in temperature for T e > 1.5 keV gas (Mohr et al. 1999). In theory, 
cluster temperature profiles may significantly affect the distance determinations through the SZE since AT cx J n e T e d£. 
The spatial filtering of the interferometer makes our SZE observations insensitive to angular scales larger than a few 
arcminutcs. Therefore we are relatively insensitive to large scale temperature gradients. However, we are sensitive to 
temperature gradients on smaller scales, for example at the center of cooling flow clusters. 

A mixture of simulations and studies of nearby clusters suggests a 10% effect on the Hubble parameter (Inagaki et al. 
1995; Roettiger et al. 1997) due to departures from isothermality. The spatial filtering of the interferometer is not accounted 
for in these studies and thus provides a conservative estimate. We include a conservative ±10% effect on the inferred 
Hubble parameter due to departures of isothermality, consistent with both cooling flow (see §6.1.3) and non-cooling flow 
departures. 

6.1.3. Cooling Flows 

Cooling flows affect the emission weighted mean temperature and enhance the X-ray central surface brightness (see, 
e.g., Fabian 1994; Nagai et al. 2000). When the cooling time at the center of the cluster is less than the age of the cluster, 
then the central gas has time to cool. This is known as a cooling flow (e.g., Fabian 1994). The cluster temperature is 
expected to decrease towards the center of the cluster, which has recently been seen with both Chandra (e.g., Markevitch 
et al. 2000; Nevalainen et al. 2000) and XMM-Newton (e.g., Tamura et al. 2001). 

A characteristic cooling time for the ICM is the available radiative energy divided by its luminosity given by 

3kT e n to t _ 3fcT e fi a 
tcooi ~ 777 — 777 j (,yj 

2An e n H 2An e [i tot 

where A is the bolometric cooling function of the cluster and all quantities are evaluated at the center of the cluster. 
Cooling flows may occur if the cooling time is less than the age of the cluster, which we conservatively estimate to be the 
age of the universe at the redshift of observation, t coo i < tjj{z). 

As a check, we calculate t coo i/tjj ratios for each cluster analyzed by Mohr et al. (1999). We check our cooling flow 
and non-cooling flow determinations versus those of Peres et al. (1998) and Fabian (1994). Of the 45 clusters in the 
Mohr sample, 41 have published mass deposition rates. We assume the cluster does not contain a cooling flow if its mass 
deposition rate is consistent with zero, otherwise it is designated as a cooling flow cluster. We are able to predict whether 
a cluster has a cooling flow or not with a 90% success rate, suggesting that the ratio t coo i /tn presented in equation (9) is 
a good predictor for the presence of a cooling flow. 

The ratio t coo iltn{£) for each cluster is summarized in Table 10 for the same three cosmologies used to determine the 
Hubble constant. The central densities, n e o, used in this calculation are determined by eliminating D A in equations (2) 
and (4) in favor of n e0 . From this analysis, the seven clusters R1347, MS1358, A611, A1835, A2261, A1689, and A1413 are 
cooling flow clusters. Both A1995 and A1413 are borderline cases. Such clusters are expected to have falling temperatures 
towards the center of the cluster. 

In principle, the multiphase medium expected in cooling flow clusters could introduce large biases in isothermal beta 
model SZE and X-ray distances (Nagai & Mohr 2002). However, recent observations with Chandra and XMM-Newton 
suggest that cooling flows are not as strong as previously expected (e.g., Fabian et al. 2001; Peterson et al. 2001). To 
estimate the possible effects of cooling flow-like temperature profiles, we adopt the deprojected temperature profile of A1835 
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from an analysis of XMM-Newton observations and determine the change in the angular diameter distance introduced by 
this profile. The inclusion of the A1835 temperature profile reduces the angular diameter distance by ~ 10%, causing a 
~ 10% underestimate in the Hubble constant from cooling flow clusters when an isothermal analysis is performed (A. D. 
Miller 2001, private communication). In addition, a theoretical examination of the effects of cooling flows on SZE and 
X-ray determined distances suggests a ~ 10% underestimate of H from an isothermal analysis (Majumdar & Nath 2000). 
Assuming all seven cooling flow clusters in our sample produce a similar 10% bias in Hq, the average underestimate in 
H for our 18 cluster sample is ~ 4%. 

We combine the uncertainty from cooling flow and non-cooling flow (see §6.1.2) departures from isothermality into a 
conservative ±10% effect on the inferred Hubble parameter. 

6.1.4. Clumping; Small Scale Structure 

Clumping of the intraclustcr gas is a potentially serious source of systematic error in the determination of the Hubble 
constant. Unresolved clumps in an isothermal intracluster plasma will enhance the X-ray emission by the factor 

C^M.. (10) 

( n e) 

The cluster generates more X-ray emission than expected from a uniform ICM leading to an underestimate of the angular 
diameter distance (D A oc 5" 1 ) and therefore an overestimate of the Hubble parameter for C > 1. Unlike the orientation 
bias which averages down for a large sample of clusters, clumping must be measured in each cluster or estimated for an 
average cluster. Theoretical estimates of C are difficult because they must account for the complicated processes that 
both generate and damp density enhancements, such as preheating and gas-dynamical processes. 

There is currently no observational evidence of significant clumping in galaxy clusters. If clumping were significant and 
had large variations from cluster to cluster, we might expect larger scatter than is seen in the Hubble diagrams from SZE 
and X-ray distances (Figure 7; see also Birkinshaw 1999). Gas-dynamical cluster simulations provide an opportunity to test 
the effects of observing strategy and cluster structure on our distance determinations. These simulated clusters exhibit 
X-ray merger signatures consistent with those observed in real clusters and, presumably, they exhibit the appropriate 
complexities in their temperature structure as well. Preliminary work indicates that temperature profiles do not introduce 
a large error on our distances but that clumping in the ICM may bias distances low by up to ~ 20%. As mentioned 
above, there is no observational evidence of clumping within the ICM. However, merger signatures are common (Mohr 
et al. 1995), and the mergers are the driving mechanism behind these fluctuations in the simulated clusters (Mathiesen 
et al. 1999). 

Clumping causes an overestimate of Hq, so we include a conservative one-sided —20% possible systematic due to 
clumping. 



Table 10 
Ratio of t coo i/t H (z) 



Cosmology (il M , Ha) 



Cluster 


(0.3, 0.7) 


(0.3, 0.0) 


(1.0, 0.0) 


MS1137 


1.1 


1.9 


1.4 


MS0451 


1.0 


1.6 


1.2 


C10016 


1.5 


2.5 


1.9 


R1347 


0.1 


0.2 


0.1 


A370 


3.6 


5.8 


4.4 


MS1358 


0.4 


0.7 


0.5 


A1995 


0.9 


1.5 


1.2 


A611 


0.4 


0.7 


0.5 


A697 


1.1 


1.7 


1.3 


A1835 


0.1 


0.2 


0.2 


A2261 


0.3 


0.5 


0.4 


A773 


1.5 


2.4 


1.8 


A2163 


1.3 


2.0 


1.6 


A520 


2.1 


3.3 


2.6 


A1689 


0.3 


0.5 


0.4 


A665 


1.2 


1.8 


1.4 


A2218 


1.4 


2.2 


1.7 


A1413 


0.6 


0.9 


0.7 
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6.2. Possible SZE Contaminants 
6.2.1. Possible Undetected Point Sources in the Field 

Undetected point sources near the cluster center mask the central decrement, causing an underestimate in the magni- 
tude of the decrement and therefore an underestimate of the angular diameter distance. The synthesized beam shapes, 
which include negative sidelobes, allow both underestimates and overestimates in the magnitude of the decrement. As a 
conservative estimate of our detection threshold, we use 3 times the rms of the high reslution map, applying a > 2000 A 
cut on the baselines for each cluster data set. Placing a point source with flux equal to our detection limit near the cluster 
center and re-analyzing to find the change in the central decrement provides an estimate of the upper bound of the effects 
of undetected radio point sources. 

We have additional information on the distribution of point sources in all our cluster fields from observations at lower 
frequencies. Sources with flux densities greater than 2 mJy at 1.4 GHz appear in the NVSS catalog (Condon et al. 1998). 
We use the NVSS catalog to find point sources within 400" of each cluster center. We extrapolate the NVSS sources in 
our fields to 28.5 GHz using the average spectral index of radio sources in galaxy clusters a = 0.77 (Cooray et al. 1998), 
where 5„ oc v~ a . Extrapolated NVSS sources with fluxes greater than our 3 a threshold are ruled out by the 30 GHz 
data and their fluxes are fixed at the maximal 30 GHz 3 a value. The extrapolated NVSS sources are added to to the 
30 GHz visibilities data, which is re-analyzed, not accounting for the additional point sources. The uncertainty on the 
angular diameter distance from undetected point sources is summarized in Table 11 for each cluster field. The average 
over the 18 cluster fields yields a ~ 12% uncertainty on the Hubble parameter. 

We know that clusters have central dominant (cD) galaxies, which are often radio bright. Therefore it is likely that 
there is a radio point source near the center of each cluster. To estimate the effects of cD galaxies on the central decrement 
we pick three clusters for which we do not detect a central radio point source, A697, A2261, and A1413. We add a point 
source fixed at the optical position of the cD (Crawford et al. 1999) and vary both the flux of the cD galaxy and the 
central decrement, keeping the ICM shape parameters fixed at their best-fit values. The cD fluxes are all consistent with 
zero and the corresponding changes in the central decrement are < 2%. This suggests that undetected cD galaxies do 
not contribute significantly to the uncertainty on the Hubble constant, < 4%, within our uncertainty budget for possible 
undetected point sources. 

Table 11 

Effects of Undetected Point Sources on D a (%) 



Cluster 


Effect (%) 


MS1137 





MS0451 





CL0016 


2 


R1347 a 


32 


A370 


6 


MS1358 b 





A1995 


6 


A611 b 





A697 


(i 


A1835 





A2261 a 


46 


A773 


6 


A2163 a 


22 


A520 a 


16 


A1689 a 


18 


A665 a 


12 


A2218 





A1413 


42 


Average 


12 



a Required 30 GHz 3 a 
truncation. 

b No NVSS sources in 
the cluster field. 
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6.2.2. Kinetic SZE 

Cluster peculiar velocities with respect to the CMB introduce an additional CMB spectral distortion known as the 
kinetic SZE. The kinetic SZE is proportional to the thermal effect but has a different spectral signature so it can be 
distinguished from the thermal SZE with spectral SZE observations. For a 10 keV cluster with a line-of-sight peculiar 
velocity of 1000 km s -1 , the kinetic SZE is ~ 11% of the thermal SZE at 30 GHz. Watkins (1997) presented observational 
evidence suggesting a one-dimensional rms peculiar velocity of ~ 300 km s -1 for clusters, and recent simulations found 
similar results (Colberg et al. 2000). With a line-of-sight peculiar velocity of 300 km s _1 and a more typical 8 keV cluster, 
the kinetic SZE is ~ 4% of the thermal effect, introducing up to a ~ ±8% correction to the angular diameter distance 
computed from one cluster. When averaged over an ensemble of clusters, the effect from peculiar velocities should cancel, 
manifesting itself as an additional statistical uncertainty similar to the effects of asphericity. Therefore, we include a 2% 
(= 8 / a/T8) effect from the kinetic SZE for our 18 cluster sample. 

6.2.3. CMB Primary Anisotropics 

CMB primary anisotropics have the same spectral signature as the kinetic SZE. Recent BIMA observations provide 
limits on primary anisotropies on the scales of the observations presented here (Dawson et al. 2001; Holzapfel et al. 
2000). They place a 95% confidence upper limit to the primary CMB anisotropies of AT < 19 fiK at I ~ 5500 (~ 2' 
scales). Thus primary CMB anisotropies are an unimportant (< 2%) source of uncertainty for our observations. At 68.3% 
confidence, primary CMB anisotropics contribute a < 1% uncertainty in the measured Hubble parameter. In addition, 
CMB primary anisotropy effects on the inferred Hq should average out over the sample; with an 18 cluster sample CMB 
primary anisotropy contributes < 1% uncertainty to Hq. 

6.2.4. Radio Halos 

The SZE decrement may be masked by large scale diffuse non-thermal radio emission in clusters of galaxies, known 
as radio halos. If present, radio halos are located at the cluster centers, have sizes typical of galaxy clusters, and a 
steep radio spectrum a ~ 1 — 3 (Kempner & Sarazin 2001; Giovannini et al. 1999; Moffet & Birkinshaw 1989; Hanisch 
1982). Similar structures at the cluster periphery, usually with an irregular shape, are called relics. In general, radio 
halos and relics are rare phenomena that are present in rich, massive clusters, characterized by high X-ray luminosity and 
high temperature (Giovannini & Feretti 2000; Giovannini et al. 1999). Cooling flow clusters rarely contain radio halos. 
Because halos and relics are rare, little is known about their nature and origin but they are thought to be produced by 
synchrotron emission from an accelerated or reaccelerated population of relativistic electrons (e.g., Jaffe 1977; Dcnnison 
1980; Roland 1981; Schlickeiser et al. 1987). Shocks from cluster mergers may be the acceleration mechanism though 
there are numerous theories (e.g., Jaffe 1977, 1980; Dcnnison 1980; Roland 1981; Schlickeiser et al. 1987; Ensslin et al. 
1998; Blasi & Colafrancesco 1999; Dolag & Ensslin 2000; Liang et al. 2000). 

According to the literature (Kempner & Sarazin 2001; Giovannini & Feretti 2000; Giovannini et al. 1999), the following 
clusters in our sample exhibit radio halos: C10016, A773, A2163, A520, A665, A2218. In Figure 9, we show NVSS 1.4 GHz 
image contours (Condon et al. 1998) overlaid on color scale images of the SZE cluster emission. Contours are multiples of 
twice the rms in the NVSS maps (rms ~ 0.45 mJy beam -1 ). It is apparent that many of these halos are at the 2 a « 0.9 
mJy beam -1 level. The brightest known halo is seen in A2163 with a peak brightness of ~ 5.4 mJy beam -1 . 

For each of the halo clusters, we conservatively model the halo as a point source at the cluster center with flux from 
an a = 1 extrapolation of the peak NVSS halo flux and re-analyze to determine the effect on the central decrement. The 
average effect on the central decrement from the clusters with radio halos in our sample is 4%, excluding A2163 which 
shows a ~ 10% effect on the central decrement. Radio halos typically have spectral indices a > 1.5, making the a = 1 
extrapolation a conservative upper bound for the effects of radio halos. Averaged over the entire 18 cluster sample, the 
a = 1 extrapolation results imply a ~ 3% overestimate (—3% effect) on our inferred Hubble parameter (Hq cx AT -2 ) 
from radio halos. 

6.2.5. Imprecisely Measured Primary Beam 

The primary beam is determined from holography measurements at both OVRO and BIMA. The effect of the primary 
beam on the interferometric observations is a convolution in the Fourier plane (see §3.1.1) or equivalently, an attenuation 
across the field of view. Therefore, differences in primary beam shape simply alter the smoothing kernel in the u-v plane 
only slightly, having a small affect on the derived distances. 

To assess the effects of the primary beam quantitatively, we fit to an OVRO and BIMA data set using a cluster model 
attenuated by a Gaussian beam with different assumed FWHMs. There is no significant change in the central decrements 
when using a Gaussian approximation for the primary beam instead of the measured beam. Even with an unrcalistically 
large ±0.2' uncertainty in the primary beam FWHM, the uncertainty introduced in the Hubble constant is < 3% (< 2% 
in ATo). Artificially broadening the wings of the real beam has a negligible effect on the derived central decrements. We 
adopt a 3% uncertainty in Hq as a conservative estimate of the maximum effects from an imprecisely measured primary 
beam. 
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Fig. 9. — -SZE from our own 30 GHz observations (color scale) with NVSS 1.4 GHz contours. The 30 GHz images are 
the color scale version of the SZE shown in Figure 2. The color scale wedge above each image shows the range in the 
flux density of the 30 GHz map in units of mJy beam -1 . Contours arc multiples of twice the rms and the NVSS rms is 
~ 0.45 mJy beam -1 . The FWHM of the 30 GHz synthesized beam is shown in the lower left hand corner of each panel 
and the 45" FWHM beam of the NVSS survey is shown in the lower right hand corner of each panel. A2163 exhibits the 
brightest 1.4 GHz radio halo in our 18 cluster sample. 
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7. DISCUSSION AND CONCLUSION 

We perform a maximum- likelihood joint fit to centimeter-wave interferomctric SZE and ROSAT X-ray (PSPC and HRI) 
data to constrain the ICM parameters for a sample of high redshift clusters of galaxies. We model the ICM as a spherical, 
isothermal (3 model. From this analysis we determine the distances to 18 galaxy clusters. Together, these distances imply 
a Hubble constant of 

( 60±|±i| km s" 1 Mpc" 1 ; O M =0.3, ft A =0.7, 
H = < 56^ t\r km s" 1 Mpc" 1 ; Q M =0.3, Cl A =0.0, (11) 
1 54 ±3 -16 km s_1 Mpc" 1 ; fi M =1.0, ft A =0.0, 

where the uncertainties are statistical followed by systematic at 68.3% confidence. The systematic uncertainties have been 
added in quadrature and include an 8% (4% in AXb) uncertainty from the absolute calibration of the SZE data, a 10% 
effective area uncertainty for the PSPC and HRI, a 5% uncertainty from the column density, a 5% (~ 20/VT8) uncertainty 
due to asphericity, a 10% effect for our assumptions of isothermality, a one-sided —20% effect from possible small-scale 
clumping in the ICM, a 12% uncertainty from undetected radio sources, a 2% (~ 8/vT8) uncertainty from the kinetic 
SZE, a 1% uncertainty from primary CMB anisotropics, a —3% effect from radio halos, and a 3% effect from an imprecisely 
measured primary beam. These systematic uncertainties are summarized in Table 9. We adopt conservative assumptions 
when gauging the effects of possible systematics. The contributions from asphericity, kinetic SZE, and primary CMB arc 
expected to average out for a large sample. 

The measured distances plotted in Figure 7 with the theoretical relation show that A370 is the largest outlier from the 
theoretical angular diameter distance relation and MS1137 has the largest distance uncertainty. A370 exhibits an almost 
2-to-l axial ratio in the knotty N-S elongation of its X-ray image (see Fig 2). An optical study of this cluster and its 
member galaxies shows that the cluster is dominated by two giant elliptical galaxies with a projected separation of about 
40", roughly in the north-south direction (Mcllicr ct al. 1988). In addition, gravitational lens models suggest that A370 
has a bimodal mass distribution with the two components separated in a roughly north-south direction (Kneib ct al. 1993; 
Soucail et al. 1988; Smail et al. 1996). The spherical model used is clearly insufficient for the complex structure of this 
cluster. The uncertainties on the distance to MS1137 are particularly large since MS1137 resides in the distant universe 
(z = 0.78), making it difficult to collect large numbers of X-ray photons. Therefore, the uncertainties on the X-ray driven 
quantities are large; in particular, the uncertainty on the measured X-ray temperature is the main contributor to the 
large uncertainty on the angular diameter distance to MS1137. The large uncertainties on the distances to both A370 
and MS1137 mean that those clusters contribute little weight in the determination of the Hubble constant. 

As discussed in §2, target clusters were originally chosen from a limited sample of known X-ray clusters. We construct 
subsamples of our cluster sample to explore the robustness of our result and to look for possible biases in our Hq 
determination. A description of the subsamples, the number of clusters in each subsample N, the Hubble constant from 
each subsample, and the x 2 an d reduced x 2 f° r the Hubble parameter are summarized in Table 12 for each subsample 
considered. Only the Qm = 0.3, = 0.7 cosmology is considered in this study. Other cosmologies will have similar 
changes to the best fit Hubble parameter. Table 12 shows that excluding the largest outlier, A370, has a negligible effect 
on the determined Hubble parameter. We also split the sample up based on having a cooling flow, by redshift, presence 
of a radio halo, based on point sources, based on right ascension and declination, based on X-ray luminosity, and based 
on membership in the EMSS survey or being an Abell cluster. 

Jones et al. (2002) constructed an orientation unbiased sample of galaxy clusters for a SZE and X-ray determination 
of the Hubble parameter. Most importantly, they drew clusters present in both the BCS (Ebeling et al. 2000a; Crawford 
et al. 1999; Ebeling et al. 1998, 1997) and NORAS (Bohringer et al. 2000) surveys with L* > 8 x 10 44 erg s" 1 and well 
above the survey flux limits (> 5 x 10~ 12 erg s" 1 cm" 2 ). Eleven clusters satisfy these criterion as well as their redshift 
range choice (0.14 < z < 0.30) and declination constraint (> 2°). They find five of the eleven to be sufficiently free from 
point sources at 15 GHz for SZE measurements. Six of these eleven arc part of our sample, which we call our orientation 
unbiased subsample in Table 12. 

All of these subsamples yield a Hubble constant with 1 a statistical uncertainties consistent with the Hq from the entire 
18 cluster sample. This argues in favor of a robust Hq determination. 

We compare our results with other SZE determined distances to clusters in our sample in Table 13. Only statistical 
uncertainties are included. There are nine clusters in our sample that also have previously determined SZE distances. All 
of the 1 a confidence regions agree with our own, with the exception of C10016, R1347, A773, and A665. The systematic 
uncertainties on the angular diameter distance are > 30% for one galaxy cluster. Therefore, all of the distances are in 
reasonably good agreement, even after accounting for shared systematics (namely most use ROSAT X-ray data and ASCA 
X-ray temperatures). 

Many of the systematics can be approached and reduced through improved observations. For example, Chandra and 
XMM-Newton are now producing temperature profiles of galaxy clusters (e.g., Markevitch et al. 2000; Ncvalainen et al. 
2000; Tamura et al. 2001). The unprecedented angular resolution of Chandra will provide insight into possible clumping 
in clusters. The effects of undetected point sources are being addressed with multi-wavelength (5 and 8 GHz) VLA 
observations of many of our cluster fields. In addition, there is a project to produce a ~ 1% calibration of the modified 
OVRO and BIMA SZE systems and the current generation of X-ray satellites will reduce the X-ray absolute calibration 
uncertainty to the few percent level. 

The 18 cluster distances presented here are beginning to probe the shape of the angular diameter distance relation. 
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Table 12 

Ho from Different Subsamples for the A Cosmology 







H 






Subsamplc 


N 


(km s" 1 Mpc" 1 ) 


x 2 


2 a 
A- red 


All 


18 


60l 4 


16.5 


0.97 


No A370 


17 


6111 


10.8 


0.68 


No MSI 137 or A370 


16 


Qltl 


9.8 


0.65 


No cooling flow 


11 


59 + t 


14.5 


1.45 


Only cooling flow 


7 




1.7 


0.28 


z > 0.27 


9 


62^ 


9.7 


1.21 


z < 0.27 


9 


59 +6 5 


6.6 


0.83 


No point sources 


4 


53t? 


3.4 


1.13 


Only point sources 


14 


62+| 


12.2 


0.94 


Only radio halo 


6 




6.7 


1.34 


No halo 


12 


62 + s 


9.3 


0.85 


RA< 11. 5 h 


9 


61±S 


13.6 


1.70 


RA> 11. 5 h 


9 


60+° 


2.9 


0.36 


Deo 30° 


9 




8.3 


1.04 


Dec< 30° 


9 


6i±i 


8.1 


1.01 


Lk > 14.5 x 10 44 erg s- 1 


9 


63 + s 


4.4 


0.55 


Zx < 14.5 x 10 44 erg s" 1 


9 


57+? 


11.5 


1.44 


EMSS clusters 


5 


66+£ u 


3.5 


0.88 


Abcll clusters 


13 


581^ 


12.2 


1.02 


Orientation unbiased 


6 


58±» 


6.2 


1.24 



a Xred = X 2 /(N — 1), where N — 1 is the degrees of freedom. 
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Table 13 

Comparison of SZE Determined Distances 









Cluster 


(Mpc) 


reference 


C10016 


2041 ± 499 


This work 




1788 ± 664 


H98 




1100 ±295 


G02 


R1347 


1221 ±356 


This work 




1890 ± 644 


P01 




1897 ±401 


K99 


Aoy I 


yys ± z < 4 


liiis work 




1044 ± 239 


J02 


A1835 


1027 ± 196 


This work 




867 ±411 


MOO 


A773 


1450 ± 347 


This work 




1002 ±257 


S99 


A2163 


828 ± 193 


This work 




728 ± 387 


H97 




615 ±327 


L98 


A665 


466 ± 198 


This work 




1017 ±229 


B91 


A2218 


1029 ± 346 


This work 




616 ± 118 


B94 




720 ± 422 


T98 




1201 ±343 


J02 


A1413 


573 ± 161 


This work 




565 ± 164 


G02b 



a Includes approximate 68.3% con- 
fidence statistical uncertainties only. 

REF: B91-Birkinshawet al. 1991; 
B94-Birkinshaw & Hughes 1994; 
G02-Grainge et al. 2002a; G02b- 
Grainge et al. 2002b H97-Holzapfel 
et al. 1997b; H98-Hughes & Birkin- 
shaw 1998; J02- Jones et al. 2002; 
K99-Komatsu et al. 1999; L98- 
Lamarre et al. 1998; MOO-Mauskopf 
et al. 2000; POl-Pointecouteau et al. 
2001; S99-Saunders et al. 1999; T98- 
Tsuboi et al. 1998; 
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Moreover, constructing subsamples from our 18 cluster sample based on such considerations as cooling flows, rcdshift, and 
X-ray luminosity, does not significantly affect the best fit Hq, suggesting a robust determination of the Hubble parameter. 
Systematics currently dominate the uncertainty in our determination of the Hubble parameter. These systematics can 
and will be addressed with current radio observatories (OVRO, BIMA, and VLA) and X-ray satellites (Chandra and 
XMM-Newton) . With a sample of high redshift galaxy clusters, this method can be used to constrain the geometry of the 
universe, providing a valuable independent check of the recent supernovae type la (Schmidt et al. 1998; Riess et al. 1998; 
Perlmutter et al. 1999; Goobar et al. 2000) and primary CMB power spectrum results (Pryke et al. 2002; de Bernardis et al. 
2002; Stompor et al. 2001). We emphasize that SZE and X-ray determined distances are independent of the extragalactic 
distance latter and do not rely on clusters being standard candles or rulers. 

A complete review of other distance determination methods is beyond the scope of this paper. We will just touch 
on a few methods that complement SZE determined distances. The SZE derived distances are direct, making them an 
interesting check of the cosmological distance ladder. Our measurement of Hq in the distant universe, agrees within the 
uncertainties with the Hubble Space Telescope (HST) Hq Key Project results, which probes the nearby universe. The 
HST Hq Key Project finds Hq = 72 ± 3 ± 7 km s _1 Mpc -1 (Freedman et al. 2001), where the uncertainties are statistical 
followed by systematic at 68.3% confidence. Though few in number, there are other methods that yield distances that are 
independent of the extragalactic distance ladder. Recent observations of masers orbiting the nucleus of the nearby galaxy 
NGC 4258 (Herrnstein et al. 1999) illustrate a method of determining direct distances in the nearby universe. Time delays 
produced by lensing of QSO's by galaxies are another direct distance indicator that can probe the high- redshift universe 
(for recent examples, see Fassnacht et al. 1999; Biggs et al. 1999; Lovell et al. 1998; Barkana 1997; Schechter et al. 1997). 

SZE surveys provide a promising method of detecting high redshift galaxy clusters (e.g., Holder et al. 2000; Barbosa 
et al. 1996). These surveys will provide large catalogs of high rcdshift galaxy clusters required to determine the geometry 
of the universe from SZE and X-ray determined direct distances. 

This work is the thesis work of E. D. R and benefited from useful discussions from very many people. In particular, 
we thank Gilbert Holder, Jack Hughes, Carlo Graziani, Sandeep Patel, and Thomas Crawford for useful and stimulating 
discussions. This project would not have been possible without all of the help and support of both the OVRO and BIMA 
staff over many years. In particular, E. D. R would like to thank D. Woody, S. Scott, R. Lawrence, J. R. Forster (Rick), 
R. Plambeck (Dick), J. Wirth (Red), M. Warnock, and M. Masters for both helping to keep our observations going and 
for teaching E. D. R. all about the details of each observatory. Dave, Rick, and Dick have been especially patient teachers. 
We would also like to thank the staff of the High Energy Astrophysics Science Archive Research Center (HEASARC). 
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